Abstract 


Presentation 


Table  of  Contents 


Next  Abstract 


Simulation  and  Analysis  of  Adaptive 
Interference  Suppression  for  Bistatic 
Surveillance  Radars* 


Fred  Pearson 
Geordi  Borsari 

13  March  2001 


*This  work  was  sponsored  under  United  States  Air  Force  Contract  F 1 9628-00-C-002.  Opinions,  interpretations,  conclusions 
and  recommendations  are  those  of  the  authors  and  are  not  necessarily  endorsed  by  the  United  States  Air  Force. 

MIT  Lincoln  Laboratory 

999999-1 
CFP  09  Jan  2001 


Report  Documentation  Page 

Form  Approved 

0MB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 

1.  REPORT  DATE 

14  MAR  2001 

2.  REPORT  TYPE 

N/A 

3.  DATES  COVERED 

4.  TITLE  AND  SUBTITLE 

Simulation  and  Analysis  of  Adaptive  Interference  Suppression  for 

Bistatic  Surveillance  Radars 

5a.  CONTRACT  NUMBER 

F19628-00-C-0002 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

Fred  Pearson;  Geordi  Borsari 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

MIT  Lincoln  Laboratory  244  Wood  St  Lexington,  MA  02420 

8.  PEREORMING  ORGANIZATION  REPORT 

NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

See  ADM001263  for  entire  Adaptive  Sensor  Array  Processing  Workshop. 

14.  ABSTRACT 

See  briefing  charts 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIEICATION  OE: 

17.  LIMITATION 

OE  ABSTRACT 

uu 

18.  NUMBER 
OE  PAGES 

21 

19a.  NAME  OE 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


Outline 


Problem  Overview 


*  Bistatic  Algorithms  -  Description  and  Analysis 


*  Summary  and  Future  Work 
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Problem  Overview 


*  Bistatic  geometry  invoives  separate  transmit  and  receive  piatforms 

-  Platforms  are  moving  independently 

*  Receive  oniy  piatform  for  surveiiiance  or  strike 

-  Extend  coverage  area 

-  Improve  target  localization 

-  No  transmitter  on  receive  platform 

Reduce  size,  weight,  power 
Improve  stealthiness 
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Challenges  for  Bistatic  Operation 
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*  Benefits  of  bistatic  operation  come  at  a  price 

-  Azimuth  /  Doppler  structure  of  clutter  interference  varies  with  range 

*  Chaiienge  is  to  find  training  strategies  to  estimate  covariance  R 
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Algorithm  Development  Approach 
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Test  ideas  in  idealized  geometry 
(Covariance  analysis) 


Identify  modifications  as  needed 


Identify  reduced  number  of 
promising  approaches 


Test  for  robustness  in  more 
realistic  situations 

(Training  strategies  /  simulated  data ) 


*  Covariance  model  is  used  to  compare  algorithms  with 

-  large  number  of  geometries 

-  coarse  range  sampling 

*  Modeling  goal  is  to  quickly  survey  algorithm  performance 

-  simplified  scattering  model 

*  Time  series  model  is  used  to  compare  algorithms  with 

-  small  number  of  geometries 

-  fine  scale  range  sampling 

*  Designed  to  examine  “real  world”  effects  on  algorithm  performance 
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Outline 


*  Problem  Overview 

*  Bistatic  Algorithms  -  Description  and  Anaiysis 

^  Algorithm  description 

-  “Standard”  2  -  bin  Post  -  Doppler 

-  2  -  bin  Post  -  Doppler  with  Derivative  Based  Updating  (DBU) 

-  Uses  only  radar  data  but  doubles  the  degrees  of  freedom  (DOF’s) 

-  Requires  increased  sample  support 

-  2  -  bin  Post  -  Doppler  with  High  Order  Doppler  Warping  (HODW) 

-  Uses  knowledge  of  bistatic  ciutter  ridge 

-  Receiver  must  know  position  and  velocity  of  transmitter 

-  Algorithm  performance 

*  Summary  and  Future  Work 
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2  -  Bin  Post  -  Doppler  Algorithm 


V 


*  Two-Bin  nulling  algorithm: 

-  Train  on  clutter  in  Doppler  bin  #’s  1  and  2  to  null  clutter  at  the  target  Doppler 
frequency 

*  Well  established  approach  for  monostatic  STAR  applications 

-  Typically  assume  range  invariance  and  estimate  covariance  with  range  average 
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Derivative-Based  Updating  Aigorithm 


*  Derivative-Base  Updating  Algorithm  (DBU): 

-  Hayward  (1996),  Zatman  &  Kogon  (2000  ASAP),  Zatman  (2001  ASAP) 


◄ - Training  Set - ► 


/  X  dw 

◄ - w(r)  =  Wn  +  r  —  ► 

V  dr 


w(r^  ) 


w(r2) 


-► 

Range 

(Time) 


*  Assumes  weight  vector  varies  iineariy  with  range 

-  Effectiveness  depends  on  accuracy  of  weight  vector  model 

*  Doubies  the  number  of  degrees  of  freedom  (DOF)  in  the  STAR 
problem 

-  Covariance  matrix  size  is  doubled 

-  Number  of  training  samples  required  to  estimate  covariance  is 
doubled 
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Derivative  Based  Updating  -  Interpretation 


•  Assume  optimal  filter  =  Wg  +  k  w’  (at  relative  range  gate) 

•  Xi^  +  kw’^^x,^  =[wj^  w’^^][Xi^;kxJ 

•  Form  sample  set  based  on  extended  vector  [x,^ ;  kx,^]  to  obtain 
extended  covariance 


Res,  =  (1/N) 
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[Wgf^  w’f^][x^;kxj=  [v»  0](R^JHx^;kx^]  =  v»  x^ 
(have  used  sample  set  symmetry  ( 51*  =  0)  and  Ri^  =  <  X/^  x^^  >  =  Rg  +  kR’) 

k 

DBU equivalent  to  applying  filter  v 

with  =  (l-k  Rg^  R’ )  (Rg -a  R’  Rg^  R’ 

•  First  order  perturbation:  '  =  ( Rg  +  k R’ ~  (I -k Rg-''  R’ )  Rg-' 

-  DBU  matches  perturbation  up  to  terms  quadratic  in  R’ 

-  the  a  term  grows  quadratically  with  the  size  of  training  set 

MIT  Lincoln  Laboratory 


999999-9 
CFP  09  Jan  2001 


Doppler 


High  Order  Doppler  Warping  (HODW) 


Clutter  Ridge  (near  target  azimuth) 
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*  In  each  Doppler  filter  apply  a  range-dependent  Doppler  frequency  shift 

-  Shift  is  different  in  each  Doppler  filter,  at  each  range 

Original  warping  algorithm  used  same  shift  in  each  Doppler  filter 

-  Interference  structure  nearly  homogeneous  In  range  for  each  output  Doppler  bln 

Clutter  ridge  calculation  requires  knowledge  of  transmitter  position  and  velocity 
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Doppler  Frequency  -  Hz 


High-Order  Doppler  Warping 

Bistatic  Space  to  Air  Example 


Clutter  Ridges  Over  6  km  at  Target  Range 
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*  Frequency  shift  is  derived  from  the  clutter  ridge  geometry 

-  Clutter  ridge  multiplicity  (front  lobe  /  back  lobe,  aliasing)  resolved  by 
choosing  highest  transmit  power  branch 

*  “High  Order”  Warping  has  made  the  clutter  interference  range 
invariant”  on  a  bin  by  bin  basis 
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Outline 


•  Problem  Overview 

•  Bistatic  Algorithms  -  Description  and  Anaiysis 

-  Algorithm  description 

-  “Standard”  2-bin  Post  -  Doppler 

-  Derivative  Based  Updating 

-  High  Order  Doppler  Warping 

^  Algorithm  performance 

•  Summary  and  Future  Work 
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Measuring  Performance 


•  Standard  measure  of  performance  is  S/A//?  Loss 

•  For  signal  element  response  vector  v  (j  v  p  =  1 )  and  filter  w : 

-  SINR  =  lsplv^vp/(v^Rw) 

where  R\s  the  true  “interference  +  noise”  covariance  matrix  <  xx^  > 
and  s  is  the  signal  amplitude 

•  For  uncorrelated  noise  (unit  power)  <  n  n^>  =  I  and  with  w  =  v 

-  SNR  =  lsp  lv»  vp /(v»  v)  =  lsp 

•  For  correlated  noise  <  n  n^>  =  A/and  with  w  =  N-^  v 

-  SNR  =  Isp  N-^  Vp/(V»  N-^  v)  =  I  s  P  v^N'^v 

•  RsWo\sSINRLoss  =  lwHvp/((wH  Rw)(v»N V))  <1 

-  Optimal  w  =  R^  vand  max(  SINR  Loss  )  =  v/(v^N  v) 

-  In  practice  use  estimated  R^^,  and  iv  =  R^^p  v 
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Latitude  -  deg. 


Case  1 :  Air  to  Air  Geometry 
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Training  Region  -  20  km 
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*  Moderate  variation  of 
clutter  ridge  with  range 
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Range  Rate  -  m/s 


Algorithm  Performance  -  Bistatic  Air  to  Air 

(Case  1 ) 


*  Standard  Sample  Covariance  Matrix  approach  significantiy  degraded 

-  Only  moderate  variation  of  clutter  interference  structure  across  training  region 

-  Standard  approach  preserves  60%  of  useable  Doppler  space  (UDSF) 

*  Both  DBU  and  HODW  methods  yield  near  -  ideal  performance 

-  DBU  preserves  80%  UDSF,  HODW  85%,  Ideal  85% 
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Latitude  -  deg. 


Case  2:  Space  to  Air  Geometry 
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*  Clutter  ridge  varies  rapidly 
with  range 
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Range  Rate  -  m/s 


dl 


Algorithm  Performance  -  Bistatic  Space  to  Air 

Case  2 


*  Standard  Sample  Covariance  Matrix  approach  performs  poorly 

-  Very  rapid  variation  of  clutter  interference  structure  across  training  region 

-  Much  worse  performance  than  in  air  to  air  case 

-  UDSF  degrades  from  45%  with  4  km  training  to  25%  with  6km  training 

*  Both  DBU  and  HODW  methods  again  yield  near  -  ideal  performance 

-  UDSF  is  80%  for  both  DBU  and  HODW,  UDSF  for  ideal  is  90% 
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Bistatic  STAR  Aigorithms  -  Recap 


*  standard  training  approach  for  STAP  works  poorly 

-  Poor  choice  for  non  -  stationary  interference 

*  DBU  approach 

-  Advantages 

No  knowledge  of  transmitter  position  and  velocity  required 

-  Disadvantages 

Doubles  the  STAP  degrees  of  freedom 

Doubles  the  number  of  training  sampies  required 
Increases  cost  of  weight  computation  by  factor  of  8 

No  significant  impact  on  weight  application  computation 

*  HODW  Approach 

-  Advantages 

No  increase  in  degrees  of  freedom  required 
Fully  adaptive  in  spatial  dimension 

-  Disadvantages 

Requires  knowledge  of  transmitter  position  and  velocity 
Increased  complexity  of  Doppler  filtering 
FFT  techniques  may  not  be  possible 
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•  Problem  Overview 


*  Bistatic  Algorithms  -  Description  and  Anaiysis 


^Summary  and  Future  Work 
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Summary 


*  Bistatic  clutter  interference  suppression  poses  new  challenges 

-  Clutter  interference  exhibits  strongly  range  dependent  structure 

*  Doppler  warping  technique  generalized 

-  “High  Order  Doppler  Warping”  algorithm 

*  2-bin  Post-  Doppler  Algorithms  examined  both  with  covariance 
analysis  and  more  realistic  direct  time  series  analysis 

*  Preliminary  assessments  of  selected  algorithms  in  Air  -  to  -  Air 
and  Space  -  to  -  Air  bistatic  scenarios  presented 

All  algorithms  rely  on  sample  average  over  range  to  estimate  clutter 
Interference  covariance 

-  Standard  training  -  POOR 

(no  attempt  to  address  range  variation) 

-  Derivative  Based  Updating  (DBU)  -  GOOD 

Requires  doubiing  probiem  dimensionaiity 

-  High  Order  Doppler  Warping  (HODW)  -  GOOD 

Requires  knowiedge  of  transmitter  position  and  veiocity 
Doppier  fiiter  impiementation  more  compiex 
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Future  Directions 


Next  Abstract 


*  Extend  analyses  to  other  engagement  geometries 

*  Assess  impact  of  imperfections 

-  Array  element  calibration  uncertainties 

Both  DBU  and  HODW  are  fully  data  adaptive  in  the  spatial  dimension 
No  deterministic  spatial  transformations 
Anticipate  impact  simiiar  to  that  on  monostatic  STAR 

-  Engagement  geometry  uncertainties 

HODW  requires  a  priori  knowledge  of  transmitter  position  and  velocity 

*  Develop  computational  complexity  estimates  for  HODW 

-  Determine  optimal  implementation  strategy 
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